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quartz-crystal whereas each satellite clock is based upon a more accurate atomic reference 
clock. The receiver has the orbital patterns of the satellites stored in a memory, and it can 
determine the orbital position of the satellites based on the time of its clock. The receiver reads 
the timing marks on the satellite's signal, and compares them against it own clock to determine 
the transmission time from satellite to receiver. 

[0035] The satellite's low-frequency (e.g., 50Hz) information signal provides the 
least precise timing information, the C/A-code signal provides the next most precise timing 
information, and the P-code signal provides the most precise timing information. The pseudo- 
range is determined from the low-frequency information signal and the C/A-code signal for 
civilian users and some military users, and is determined from the low-frequency information 
signal and the P-code signal for most military users. Accurate use of the P-code signal requires 
knowledge of a certain code signal that is only known to military users. Better precision than 
that provided by the P-code signal can be obtained by measuring the phase of the satellite 
carrier signal in a differential navigation DN mode using two receivers. 

[0036] Fig. 1 shows a schematic drawing of an exemplary differential navigation 
system incorporating certain features of the present invention. The differential navigation 
system includes a Rover GPS unit 1 10 and a Base GPS Unit 1 14. The location of the Base 1 14 
is known and the system is attempting to determine the location of the Rover 110. As is well 
known in the art, Rover 1 10 receives signals from satellites 102, 104, 106, 108 via antenna 1 12. 
Similarly, Base 1 14 receives signals from satellites 102, 104, 106, 108 via antenna 116. The 
Base 1 14 transmits differential correction signals to the Rover 1 10 via communication channel 
118. Communication channel 1 1 8 may be, for example, a wired or wireless connection. The 
techniques for determining the location of the Rover 1 10 in the illustrative system shown in 
Fig. 1 are well known in the art, for example as described in, Bradford W. Parkinson and James 
J. Spilker Jr., Global Positioning Theory and Applications, Volume 163 of Progress In 
Astronautics and Aeronautics, published by the American Institute of Aeronautics and 
Astronautics, Inc, Washington D.C., 1996. 

[0037] A laser unit 120 associated with Base 1 14 transmits laser beam(s) 122 to an 
optical sensor 124 on the Rover 1 10 in order to precisely measure an elevation angle 126 
between reference points on the Base and Rover. A local reference plane such as 128 is used 
for measuring the elevation angle. Additional details regarding the generation, transmission, 
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fact that the corrected Rover position belongs to the surface of the cone, may be involved 
into the floating and fixed ambiguity estimation schemes. 

[0064] Let us divide the matrix D k and the vector R k in the equation (1) into blocks 
in accordance with division of the vector of unknowns Y k into three parts: 



A = 



A*,* : Ac/,* : Ac/?,* 



A,,* ; a,, 



D. 



tn t k 



\^Aw,A: 



A/,* j A*,* ) 



t,k 



\ R n,k ) 



(3) 



We will omit the epoch index k for the sake of brevity wherever this does not lead 
to misunderstanding. Given the vector h let C„ be the orthogonal projection onto the 
local horizon plane 



C„ =I-hti 



and 



C L = hh' 

is the matrix of orthogonal projection onto the direction h which is the orthogonal 
complement to C„. 

Provided there are no measurements errors, the identity 

||C |1 (X-r)||tan(^) = ||C 1 (^-r)|| 

must hold. 

Expanding last identity obtain that 

<p(X) = ((X - Tjq(X - T)} tm(0) -ti(X-T) = Q. 



(4) 



Let also introduce the quadratic penalty cost function for violation of the equation 



(4): 



®{X) = \cp{X)\ (5) 
Note now that solution of the equation (1) is equivalent to minimization of the 
quadratic function 

F(Y) = ftD k Y k = ^)|| 2 = \{D k Y k - R k ){DJ k - R k ) . 
The equalities (1) and (4) will be treated by least squares minimizing the 
'penalized' cost function 

AP D 

j-\<< ! 
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[0067] Particularly, both Lenstra-Lenstra-Lovasz transformation and 
X -decorrelation may be applied to the augmented matrix D k . 

[0068] The observation vector consists of all navigation parameters comprising 
between satellites and Rover and between satellites and Base pseudo-ranges and carrier phases. 
The observation vector is involved into calculation of the right-hand side vector R k in the 
equation (1) above. 

[0069] The state vector consists of all unknowns to be estimated: X,Y ? Z of the 
Rover, clock shift, ambiguity for all satellites and frequency bands. More precisely: X,Y,Z of 
the Rover, between Rover and Base clock shift, between Base and Rover difference of 
ambiguity for all satellites and frequency bands. The state vector estimated at the epoch k is 
denoted by Y k in the equation (1) above. 

[0070] The present invention does not replace or eliminate the "rounding of the float 
ambiguity vector to the nearest set of integers". Very generally speaking the term 'rounding 
off depends on the norm one uses to measure the closeness of two vectors - floating and 
integer ambiguity in our case. The integer ambiguity resolution, or rounding off, may be 

considered as a solution to the problem: minlA^ - Nil where N is the floating ambiguity 

vector, Z is the integer valued n-dimensional space (or integer lattice), N is the integer-valued 
vector of ambiguity to be found, ||| stands for the norm we use to measure the closeness 

between two vectors N and N . Using the ||| norm (or sum of absolute values of the vector 

components) for ||| results in the 'rounding-off in the usual sense as solution of the above 

minimization problem. We use another norm: weighted ||| 2 or \\\\ D 2 which results in the 

minimization problem min(N- N)'D( N - N) (MP). 

[0071] The present invention (as part of integer ambiguity resolution) replaces the 
matrix D with the augmented matrix D , where the matrix D appears in the equation (1), the 
matrix D is defined in the equation (8). Also the floating point vector N as part of the state 
vector Y k in (1) is replaced with the part of the state vector Y k in (7).] 
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